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ABSTRACT 



Context. The availability of asteroseismic constraints for a large sample of stars from the missions CoRoT and Kepler paves the way 
for various statistical studies of the seismic properties of stellar populations. 

Aims. In this paper, we evaluate the impact of rotation-induced mixing and thermohaline instability on the global asteroseismic pa- 
rameters at different stages of the stellar evolution from the Zero Age Main Sequence to theThermally Pulsating Asymptotic Giant 
Branch to distinguish stellar populations. 

Methods. We present a grid of stellar evolutionary models for four metallicities (Z = 0.0001, 0.002, 0.004, and 0.014) in the mass 
range between 0.85 to 6.0 M e . The models are computed either with standard prescriptions or including both thermohaline con- 
vection and rotation-induced mixing. For the whole grid we provide the usual stellar parameters (luminosity, effective temperature, 
lifetimes, ... ), together with the global seismic parameters, i.e. the large frequency separation and asymptotic relations, the frequency 
corresponding to the maximum oscillation power v mllx , the maximal amplitude A max , the asymptotic period spacing of g-modes, and 
different acoustic radii. 

Results. We discuss the signature of rotation-induced mixing on the global asteroseismic quantities, that can be detected observa- 
tionally. Thermohaline mixing whose effects can be identified by spectroscopic studies cannot be caracterized with the global seismic 
parameters studied here. But it is not excluded that individual mode frequencies or other well chosen asteroseismic quantities might 
help constraining this mixing. 

Key words. Asteroseismology - Instabilities - Stars: rotation - Stars: evolution - Stars: interiors 



1. Introduction 

In the recent years important efforts were devoted to improve 
our understanding of the physics of low- and intermediate- 
mass stars, and in particular to explain the abundance anoma- 
lies they exhibit along their lifetime. Rotation was shown to 
change their internal dynamics due to the transport of angular 
momentum and chemical species through the action of merid- 
ional circulation and shear turbulence, combined possibly with 
other processes induced by i nternal gravity w a ves or magnetic 
fields ( see e.g. IZahnl 1 1 992t Iz^in et al.l Il997b iMaeder & Zahnl 



1998t lTalon & Charbonnel 119981: lE ggenberger et al.l l2005t 



Charbonnel & Talon 20051 120081) . Rotation-induced mixing re- 



sults in variations of the stellar chemical properties that success- 
fully explain many o f the abundance pat t erns observed at the 
surfac e of these stars (|Palacios et a l. 2003; Charbonnel & Talonl 



120081: ISmiljanic et al.l l2010t ICharbonnel & Lagardel ^oToT 
Additionally thermohaline mixing driven by 3 He-burning was 
proposed to be the dominant process that further modifies the 
photospheric composition o f bright low-mass red giant stars 
(ICharbonnel & Zahnl I2007b1: for references on the abundance 
anoma lies at that evolution p hase see ICharbonnel & Lagardel 
(1201 Oh hereafter Paper I, and lLagarde et al.l (1201 ll) . hereafter 
Paper II). During the thermal-pulse phase on the asymptotic gi- 
ant branch (TP-AGB), thermohaline mixing was found to lead 



to lithium production (Paper I ; Stancliffel 1201 Ol) . accounting 
for Li abundances observed in oxygen-rich AGB variables of 
the Galactic disk (lUttenthaler & Lebzelterll2010l) . In summary 
and as discussed in Paper s I and II of this series (see also 
ICharbonnel & Zahn 2"007bl the effects of both rotation-induced 
mixing and thermohaline instability as described presently do 
account very nicely for most of the spectroscopic observations 
of low- and intermediate-mass stars at various metallicities. 



This has crucial consequences for the chemical evolution of 
the Galaxy (Paper II, and Lagarde et al., in preparation), and 
should also be taken into account in the other topical astrophysi- 
cal domains that use stellar models as input physics. This is par- 
ticuliarly true at the moment where asteroseismic probes bloom 
across the Hertzsprung-Russel (HR) diagram. Thanks to the re- 
cent development of dedicated satellites as Corot and Kepler, 
the internal proper t ies of stars both on the main sequence (e.g. 



tne internal properties of stars botn on tne main sequence (e.g. 
iMichel et all 120081: IChaplin et al1l2010l) and the giant branches 



(e.g. iDe Ridder et al l l2009t iBedding et all 120101) are revealed. 
Furthermore due to the high number of stars observed by these 
missions, statistical studies are possible through the determina- 
tion of global pulsation properties like the frequency of maxi- 
mu m oscillation powe r and the large freque ncy separation (see 
e.g. lMiglio et al]l2009tlChaplin et alJl201 lah . 
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In this broad context the aim of the present paper is to pro- 
vide the relevant classical stellar parameters together with the 
global asteroseismic properties of low- and intermediate-mass 
stars all along their evolution. This is done from the pre-main 
sequence (along the Hayashi track) to the early-asymptotic gi- 
ant branch (and along the TP-AGB for selected cases) for the 
grid of models computed in Papers I and II for four metallicities 
spanning the range between Z = 0.0001 and Z = 0.014 and with 
initial masses between 0.85 M and 6.0 M . This grid contains 
models computed with rotation-induced mixing and thermoha- 
line instability, along with standard models without mixing out- 
side convective regions for comparison purposes. The present 
work is a prerequisite to validate further the current theoretical 
prescriptions for the non-standard mechanisms before we test 
them through detailed seismic dissections of individual stars. 

Such a grid of stellar models at different metallicities is 
obviously a key tool for various important astrophysical top- 
ics related to e.g., stellar evolution in clusters, stellar nucle- 
osynthesis, chemical evolution, etc. They are a vailable since 
a lon g time for standard stellar models (e.g. ISchaller et al 
1992tlForestini & Charbonnelll997l:lYi et al.l2003l:ICassisi et al 



20061) . and have re cently appeared in t he literature for rotation- 



induced models dBrott et al.l 1201 ll lEkstrbm et all (12012b ). 
However these latest studies focus more on the evolution of mas- 
sive stars, and do not include thermohaline mixing nor study the 
TP-AGB phase for low- and intermediate-mass stars as do those 
we present here. 

The paper is organized as follows. In Sect. 2 we describe the 
physical inputs of the stellar evolution models. The table con- 
tent of our grids is presented in Sect. 3. Sect. 4 includes a short 
discussion on the main properties of th e models and a comp ari- 
son with the solar metallicity models of Ek strom et al.l(l2012l) . In 
Sect. 5 we present the asteroseismic parameters of models with 
and without rotation. Finally our main results are summarized in 
Sect. 6. 



2. Physical inputs 

The models are computed with the lagran gian implicit stel- 
lar evolution code STAREVOL (v3.00. See ISiess et al.1 l2000t 



iPalacios et aT1l2003l 120061: iDecressin et al.ll2009l) . In this section 
we summarize the main physical ingredients used for the present 
grid. 

2.1. Basic inputs 

The description of the stellar structure rests on the hydrostatic 
and the continuity equations, and the equations for energy con- 
servation and energy transport. To solve this system, the follow- 
ing physical ingredients are required : 

- Nuclear reaction rates are needed to follow the chemical 
changes inside burning sites, and to determine the produc- 
tion of energy by the nuclear reaction, e mw , and the energy 
loss by neutrino, e v . We follow stellar nucleosynthesis with 
a network including 185 nuclear reactions involving 54 sta- 
ble and unstable species from 'H to 37 C1. Numerical tables 
for the nucle ar reaction rates were generated from NA CRE 
compilation (lArnould et al.lll999t lAikawa et ai1l2005l) with 
the NetGen web interfacf]- 

We mainly use reaction s rates from NACRE or from 
ICaughlan & Fowled (I 19881) when NACRE rates are not avail- 



able. For proto n captures on eleme nts higher then Ne we fol- 
low ra tes from llliadis et al.l ( 120011) otherwise from lBao et al.l 
d2000h . The following reactions rates are computed by: 



3 He(D,p) 4 He ( iDescouvemont et al.ll2004T) 
3 He(aa, V) 12 C dFvnbo et al.ll2005l) 

■)4tt_ . 13MC/3 „\13r . 22 



B(/?,v)2 4 He ; 1J N0S,v) 13 C ; 22 NaG6, v) 22 Ne 
26 Alm Q3, v) 26 Mg ; 26 Alg(y6, v) 26 Mg dHoriguchi et al.l 

1996) 

14 C(p,y) 15 N JWiescheretal.il 19901) 

14 C(p,n) 14 N ( koehler&0'brierll989al) 

14 C(a,n) 17 ; 17 Q(n, 4 He) 14 C dSchatz et al.ll 19931) 
H C(a. V) 18 Q JFunck & Langankd 19891) ~ 
14 N(n,p) 14 C jKoehler & 0'brierll989bb 
14 N(p,y) 15 ( i Mukhamedzhanov et al.ll2003l) 
17 0(n,y) 18 dWagonerlll969ir~ 
22 Ne(p.T/) 23 Na (iHale et aLpX)! 
22 Ne(n,y) 23 Na (Beer et al.ll2002n 
22 Na(n,y) 23 Na ; 23 Na(a,p) 26 M, 
26 Mg (g,p) 29 Si ; 27 Al(a-,p) 30 Si 
1195 



25 



Mg(q,p) 28 Si 



jHauser & Feshbachl 



- 26 Alm(n,y) 27 Al ; 26 Alg(n,y) 27 Al JWooslev et aLlf!??! 2 ! 

The sc r eening factors are calculated with the formalism of 
iMitlerl J 19771) for weak and i nterme diate screening con- 
ditions and of iGraboske et all (1 19731) for strong screening 
conditions. 

Opacities are required to compute the radiative gradient V ra ,/ 
and the energy transport by radiative transfer. We ge nerate 
opacity tables according to llglesias & Rogersl (11996!) using 
the OPAL websitd^ for T > 8000/T that account for C and 
O enrichments. At lower temperature (T < 8000K), we use 
the ato mic and molecular opacities given by iFerguson et al.l 
(120051) . 

The equation of state relates the temperature, pressure, 
and density and thus provides different thermodynamic 
quantities (V a d, cp,. ..). In STAREVOL we follow the formal- 
ism developed b y lEggleton et al.l (1 19731) and extended by 
iPolsetalld 19951) that is based on the principle of Helmholtz 
free e nergy minimization (see lDufour|[l999l and ISiess et all 
120001 for detailed description and numerical implementa- 
tion). It accounts for the non-ideal effects due to Coulomb 
interactions and pressure ionization. 

The treatment of convection is needed to compute the 
temperature gradient inside a convective zone. It is based 
on classical mixing length formalism with omlt = 1-6, 
from solar-calibrated models without atomic diffusion nor 
rotati on computed by Geneva models (see lEkstrom et al.l 
|2012|) . We assume instantaneous convective mixing, except 
when hot-bottom burning occurs on the TP-AGB, which 
requires a ti me-dependent convective diffusi on algorithm as 
developed in lForestini & Charbonnell(ll997l) . The boundary 
between convective and radiative layers is defined with the 
Schwarzschild criterion. An overshoot parameter d ova -/H p is 
taken into account for the convective core. This parameter 
is set to 0.05 or to 0.10 respectively for stars with masses 



1 http : //www . astro . ulb . ac . be/Netgen/ form . html 



2 26 Aim and 26 Alg represent the radioactive nuclide 26 Al in its two 
isomeric states. 

3 |http : //adg . llnl . gov/Research/OPAL/opal . html | 
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below or above 2.0 MJ5 

We use a grey atmosphere where the photosphere is defined 
as the layer for which the optical depth t is between 0.005 
and 10. We define the effective temperature and radius at the 
layer where t = 2/3. 



For mass loss we use IReimersI (1 19751) formula (with t]r = 
0.5) from the ZAMS up to central helium exhaustion 



M 



LR 

-3.98 10-%— Mo.yT 1 . 



(1) 



On AGB we shift to the mass loss prescription by 
IVassiliadis & Woodl(ll993l) . 



2.2. Transport processes in radiative zones 
2.2.1. Thermohaline mixing 

Thermohaline instability develops along the Red Giant Branch 
(RGB) at the bump luminosity in low-mass stars and on the 
early-AGB in intermediate-mass stars, when the gradient of 
molecular weight becomes negative (V^ = < 0) in the 

external wing of the thin hydrogen-burning shell surro u nding 



the degenerate stellar core (ICharbonnel & Zahnl l2007blla[ LSiessI 
l2009HStancliffe et alj|2009t Icharbonnel & Lagarddl201ol) . This 



inversion of molecular weight is created by the 3 H e( 3 He, 2p) 4 He 
reaction dUlrichlll97lllEggleton et al.ll2006ll2008l) . 

The present gr id of models is computed us ing the prescrip- 
tion advocated by Charbonnel & Zahnl d2007bl) and Paper I, II. 



It is based on lUlrichl (1 19721) with an aspect ratio of instabil- 
ity fingers a - 6, in agreement with laboratory experiments 
dKrishnamurtill200'3l) . It includes the correction for non-perfect 
gas (including radiation pressure, degeneracy) in the diffusion 
coefficient for thermohaline mixing that writes: 



D, 



-v„ 



(V ad - V) 



for V,, < 0, 



(2) 



with K the thermal diffusivity ; <p = (dlnp/dlnp)pj ; 6 
-(dlnp/dlnv)pjj ; and with the non-dimensional coefficient 



n o 2 2 



(3) 



The value of a in actual stellar conditions was recently 
questionned by the results of 2D and 3D hydrodynamical 
simulations of thermohaline convection that favour a close 



to unity ( Denissenkovl 


20101 Denissenkov & Merrvfieldll201 it 


iRosenblum et al.ll201ll 


Traxler et al.ll201 ll). However these sim- 



ulations are still far from the stellar regime; therefore we decided 
to use in this Series the prescription described above since it 
successfully reproduces the abundance data for evolved stars of 
various masses and metallicities (see Papers I and II for a more 
detailed discussion). 

2.2.2. Rotation-induced mixing 

Pre-main sequence evolution along the Hayashi track is com- 
puted in a standard way (i.e., without accounting for rotation- 
induced mixing), and solid-body rotation is assumed on the 



For small-sized core its mass extent is not allowed to be larger than 



Zero Age Main Sequence (ZAMS). On the main sequence, 
the evolution of the internal angular momentum prof ile is 
accounted with the complete for malism developed by Zahnl 



dl992l) . lMaeder & Zahnl (Il998l) . and lMathis & Zahnl d2004 that 
takes into account advection by meridional circulation and 
diffusion by shear tu rbulence (see IPalacios et al.1 120031 120061: 
Decressin et al 1 120091 for a description of the implementation in 
STAREVOL). We do not take into account the inhibitory effects 
of n gradients in the treatement of rotation. 

We assume solid-body rotation in the convective regions as 
done in Paper I and Paper II, considering that the transport of an- 
gular momentum is dominated by the large turbulence in these 
regions that instantaneously flattens out the angular velocity pro- 
file as for the abundance profiles. This hypothesis leads to a min- 
imum shear mixin g approach in t he un derlyi ng radiative lay- 
ers as discussed in IPalacios et al.l d2006l) . and IBrun & Palaciosl 
d2009l) . 

On the other hand the transport of angular momentum in stellar 
radiative layers obeys an advection/diffusion equation : 



d(r 2 Q) 
dt 



= ± 2 ±^nu r) + lUs pVY f) ( 4) 

5r z or r or \ or j 



where p, r, and Q have their usual meaning. U r is the vertical 
component of meridional circulation velocity, and v v is the ver- 
tical component of the turbulence viscosity. 

The transport of chemical species resulting from meridional 
circulation and both vertical and horizontal turbulenc e is com- 
puted as a diffusive process dChabover&Zahnlll992l) . The ver- 
tical transport of a chemical species i of concentration c, is de- 
scribed by a pure diffusion equation : 



dcj 
~dt 



1 d 



+ — T \ r PAot— 
pr L or 



dci 



(5) 



diffusion processes 



where c, represents the variations of chemical composition due 
to nuclear reactions. The total diffusion coefficient D tot for chem- 
icals can be written as the sum of three coefficients: 



Ah + Aeff + D 



(6) 



with D tn the thermohaline coefficient (Sect. l2~2~TT i, D er y the effec- 
tive diffusion coeffic ient (<x 1 /£>h), and D v the vertical turbulent 
diffusion coefficient dTalon & Zahnlll997[) that is proportio nal to 
the horizontal diffusion coefficient lChabover & Zahnl 19 92). The 
expression o f the horizontal diffusion coefficient is taken from 
IZahnl(fT992l) . with its expression that prevents numerical diverge 



D h = —(\2V-aU\ 2 + U 2 ) 1/2 
C/, 



(7) 



We do not consider possible interactions between thermoha- 
line and rotation-induced mixing, nor magnetic diffusion. Under 
the present assumptions the thermohaline diffusion coefficient 
is several orders of magnitude higher than the total diffusion 
coefficient related to rotation and than magnetic diffusivity in 
the advanced phases where the thermohaline instability devel- 
ops ( see e.g., ICharbonnel & Lagarddl20iot ICantiello & Langerl 
120101 and Paper II). 

The complete treatment is aplied up to the RGB tip or up to 
the second dredge-up for star undergoing the He-flash episode. 
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2.3. Initial rotation velocity 

The initial rotation velocity of our models on the ZAMS 
is chosen at 45% of the critical velocity at that point, with 

3 1 

Vent -(f) 2 (tt) 2 • Here we take R the stellar radius computed 
without considering the stellar deformation due to rotation ; 
if we were t o take into account th e deformation of the stellar 
radius as in lEkstrom et al.l (|2012|) . then our initial velocities 
would correspond to 30% of critical velocity. This choice of 
Vmi/ Vcrit = 0.45 fits well the mean value in the observed velocity 
distribution of low- and intermediate-mass stars in young open 
clusters. This initial rotation rate leads to mean velocities on the 
main sequence between 90 and 137 km.s~'. 

We apply magnetic braking only to [1.25 M ; Z ] and 

[0.85 M r • Z = 0.0001 1 following the descriptio n of 

KawaleJd 19881) acco rding to iTalon & Charbonnell (1 19981) and 
Charbonnel & Talonl (Il999l) . 

The transport by internal gravity waves, which is efficient 
only in main sequence stars wit h effective temperatures on the 
ZAMS lower than 6500 K (see ITalon & Charbonnell 120031 is 
neglected. We do not account either for dynamo processes nor 
for the presence of fossil magnetic fields. 



2.4. Initial abundances 

Table [Upresents the initial abundances in mass fraction that we 
assume at different m etallicities. We adopt the solar mixture of 
lAsplund et aD(l2009h . except for Ne for which we use value de- 
rived bv lCunha et al.l (120061). W e use the ratio AF/AZ = 1 .29 de- 
rived by Ekstrom et al. ( 20121) to account for the enrichment in 
helium reported to enrichment in heavy elements in the Galaxy 
until the birth of the Sun. For th e primordi a l abun dances we take 
the WMAP-SBBN value from lCoc etal] d2004l) . To determine 
the initial composition of our models at a given metallicity Z we 
use the scaling (X, = X !Q ■ for all elements except for 1 Li 



Table 1. Initial abundances in mass fraction for the models at 
different metallicities 



-3 



which is taken constant (Li/H = 4.15.10" 1 ") for § < 7.8.10 
We assume [a/Fe] = at all metallicities, which has a negligi 
ble impact in the present context- 



3. Description of grids 

3.1. Content of electronic tables 

We provide files containing sets of relevant physical quantities 
as a function of time that characterize our models computed in 
the initial mass range between 0.85 M and 6.0 M with four 
metallicities Z = 0.0001, 0.002, 0.004 and 0.014 ([Fe/H]= -2.16, 
-0.86, -0.56, and respectively). For each mass and metallic- 
ity, models are computed with standard prescriptions (no mixing 
other than convection process), and with both thermohaline in- 
stability and rotational transport. For all masses the evolution is 
followed from the beginning of the pre-main sequence (along 
the Hayashi track) up to the early-AGB phase. For each model, 
we have selected 500 points to allow a good description of the 
full raw tracks. First, key evolutionary points are determined (see 
Fig.©: 

1 . beginning of the pre-main sequence; 
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06 
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-07 


4.97 
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-07 
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08 




1.80 


■ 10" 


05 
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■ 10 
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■ 10 
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1.29 
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07 


35 C1 


6.54 


■ 10" 


06 


1.87 


■ 10 


-06 
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• 10 


-07 


4.67 


10" 


08 


37 C1 


2.15 


■ 10" 


06 


6.14 


■ 10 


-07 


3.07 


■ 10 


-07 


1.53 


10" 


OS 


Others 


1.47 


■ io- 


03 


4.19 


■ 10 


-04 


2.09 


• 10 


-04 


1.05 


10" 


05 



5 Using [a/Fe] = +0.3 instead of for our [1.5 M ; Z = 0.0001] 
model leads to a decrease of the main-sequence lifetime by ~2%, and 
lowers the turnoff luminosity and effective temperature by 5% and 1% 
respectively. This is negligible compared to the effects of rotation. 



2. ZAMS defined as the time when the central hydrogen abun- 
dance X c has decreased by 0.003 in mass fraction compared 
to its initial value. 

3. turning point with the lowest T e ff on the main sequence; 

4. end of core H-burning defined as the point beyond which X c 
is smaller than 10~ 7 ; 

5. bottom of the Red Giant Branch (RGB); 

6. RGB-tip; 

7. local minimum of luminosity during central He-burning; 

8. local maximum of T e « during central He-burning; 

9. bottom of AGB : point with local minimum of luminosity 
after the loop on HB ; 

10. end of core He-burning defined as the point beyond which Y c 
is smaller than 10~ 4 . 

The data in the tables are linearly interpolated within the re- 
sults of each evolutionary model. We perform the interpolation 
as a function of time, central mass fraction of hydrogen or he- 
lium, or luminosity according to the evolutionary phase, and the 
final interpolated values are given at 499 points that we distribute 
as follows: 

- 99 points evenly distributed in time sample the pre-main se- 
quence between points 1 and 2. 
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Table 2. Description of table containing result of our evolution models, both standard case and including rotation and thermohaline 
mixing. 



Stellar parameters Surface abundances Central abundances 



- Model number 


'H 2 H 


'H 


- Maximum of temperature T max (K) 


3 He 4 He 


3 He 4 He 


- Mass coordinate of r max (M Q ) 


6 Li 7 Li 




- Effective temperature T e /r (K) 


'Be J Be 




- Surface luminosity L (L Q ) 


io B n B 




- Photospheric radius radius R^y (R Q ) 


12 C 13 C 14 C 


12 C 13 C 14 C 


- Photospheric density p ef y (g.cm -3 ) 


14 N 15 N 


14 N 


- Density at the location of T max , p max (g.cm -3 ) 


16 Q HQ 18 


16 Q n Q is 


- Stellar mass M (M e ) 


19p 


i» F 


- Mass loss rate (M Q .yr _1 ) 


20 Ne 21 Ne 22 Ne 


20 Ne 21 Ne 22 Ne 


- Age t (yr) 


23 Na 


23 Na 


- Photospheric gravity Iog(g e a ) (log(cgs)) 


24 Mg 25 Mg 26 Mg 


24 Mg 25 Mg 26 Mg 


- Central temperature T c (K) 


26 Al 27 Al 


26 Al 27 Al 


- Central pressure P c 


28 Si 


28 Si 



- Surface velocity v mr f (km.s -1 ) 

- Mass at the base of convective envelope (M ) 

- The large separation from scaling relation Av sca | e (vHz) 

- The large separation from asymptotic relation Av asymp (vHz) 

- The frequency with the maximum amplitude v max 

- The maximum amplitude A max 

- The asymptotic period spacing of g-modes All (s) 

- The total acoustic radius T (s) 

- The acoustic radius at the base of convective evelope t B cE (s) 

- The acoustic radius at the location of helium second-ionisation region t He (s) 




4.2 4 3.8 3.6 

log Teff (K) 



Fig. 1. Evolution track in the Hertzsprung-Russell diagram of 
the standard 4.0 M Q model at solar metallicity. Each phase is 
indicated by a different color : pre-main sequence (cyan), main 
sequence (black), Hertzsprung gap and red giant branch (green), 
core helium-burning phase (blue), and asymptotic giant branch 
(red). Black triangles show the points that we have selected to 
describe the tracks (see Sect J3.lt 



- 110 points evenly distributed in terms of the central hydro- 
gen mass fraction Xc sample the main sequence with 85 
points distributed between points 2 and 3, and 25 points dis- 
tributed between points 3 and 4. 



- 60 points evenly spaced in time sample the Herzsprung gap 
between points 4 and 5. 

- 80 points evenly distributed in terms of logL sample the 
RGB between points 5 and 6. 

- 20, 70 and 70 points evenly distributed in terms of Yq respec- 
tively sample the central He burning phase between points 6 
and ^| , 7 and 8 and, 8 and 9. 

- Finally, 50 points evenly spaced in logL are selected be- 
tween points 9 and 10 

For each model, we store the quantities given 
in table in a file that can be retrieved from 
http : //obswww . unige . ch/Recherche/evol/-Database- 



3.2. Comp arisons between models from \Ekstrdm et all 
120121 Geneva code) and our models (STAREVOL) 

lEkstrom et al.l (1201 2|) computed a large grid of models with ro- 
tation from 0.8 to 120 M at solar metallicity with the Geneva 
stellar evolution code. The present grid is comple mentary since 
the low-mass star models of lEkstrom et al.l (1201 2|) are computed 
only to the helium flash at the RGB tip and do not include ther- 
mohaline mixing. Also, we presently explore asterosismic diag- 
nostics during the TP-AGB. 

In order to offer this complementarity between the 2 sets of 
giants, we were careful in choosing the same input physics and 
assumptions, although some differences remain that we describe 
below. Both codes use the same physical inputs for convection 
(Schwarzschild criterion and overshoot), opacities, mass loss, 
and nuclear reaction rates, and for the mass domain explored 
the different equations of state have a negligible impact on the 
stellar structures. Besides, the initial abundances of our models 

6 For low-mass stars below 2.0 M G in which He-flash occurs we do 
not include table points between evolutionary points 6 and 7. 
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Z = 0.0001 




4.4 4.2 4 3.8 3.6 
log Teff (x) 



4.4 4.2 4 3.8 3.6 
log Teff (x) 



Fig. 2. Theoretical evolution tracks in the HR diagram for the en- 
semble of calculated models for all metallicities (Z , Z=0.004, 
Z=0.002, and Z=0.0001, from top to bottom). On left panels the 
tracks are shown for the rotating models from the pre-main se- 
quence on. On the right panels, both standard and rotating pre- 
dictions are shown (dashed and solid lines respectively) from the 
ZAMS and beyond. 



are similar apart from the larger number of species followed in 
STAREVOL. STAREVOL has indeed a more extended network 
of nuclear reactions than the Geneva code, which allows to 
follow in particular the evolution of unstable elements such as 
13 N, 15 O, and 1 Be. As a consequence, the convective cores are 
smaller on the main sequence in STAREVOL models (by ~7% 
for the [4.0 M Q , Z ] model), and the tracks in the HR-diagram 
are slightly less luminous (Ljo is ~10% lower for the [4.0 M , 
Z Q ] model). This results in a difference in the lifetime on the 
main sequence 4% for the [4.0 M , Z ] model. This is the main 



difference we could identify when comparing standard models 
computed with the two codes. 

For the rotational transport both codes follow the advection 
by meridional circulation and the diffusion by shear turbulence. 
We carefully checked that the use of a different prescription for 
the turbulent diffusion coefficient D v dTalon & Zahnl d 19971) in 
STAREVOL and lMaeden d 19971) in Geneva code) has no impact 
in the mass domain we explore as the mixing above the con- 
vecti ve core is dom inated by D e ff. Both codes use prescriptions 
from IZahnl d 19921) for the horizontal diffusion coefficient, D/,. 
However in STAREVOL we always use Eq.7 that is given by 
IZahnl d 19921) to avoid the divergence of effective diffusion co- 
efficient (oc -gj), while in Geneva the original Eq.2.29 of IZahnl 
d 19921) is used. Surprisingly Eq.7 implies higher values of Z?/, 
at the edge of the convective core (by about a factor of 4 in 
the [4.0 M , Z ] model on the main sequence), which leads to 
less mixing in these central regions. Combined with the fact that 
STAREVOL models have a lower initial rotation velocity on the 
zero age main sequence, this results in a lower lifetime and lower 
luminosity (the maximum difference being of 15% for the [4.0 
M , Z ] model). 



4. Global properties of the models 

The effects of rotation-induced mixing and thermohaline insta- 
bility on the surface chemical properties of the grid stars were 
extensively discussed in Papers I and II. Here we focus on the 
global properties of the models 



4. 1. Hertzsprung-Russell diagrams and logTc vs logp c 

The evolutionary tracks in the HR diagram are shown in Figj2] 
for all the models computed in this study, with and without ro- 
tation. We also present the evolution of the central temperature 
and density in Figf3]for the solar-metallicity case. 



4.1.1. Metallicity effects 

Let us first recall the impact of metallicity as already reported 



in the literature (e.g. ISchaller et al.lll992t iSchaereret al.lll993t 

ar 



ICharbonnell[l994t iHeger & Langerll2000t lMaederil2009l) . For a 



given stellar mass one notes the following effects when metallic- 
ity decreases, due to the lowering of the radiative opacity: 

- The ZAMS is shifted to the blue, and both the stellar lu- 
minosity and effective temperature are higher at that phase 
compared to solar metalicity models. 

- The positions of the RGB and of the AGB are shifted to the 
blue, and the ignition of helium (in degenerate conditions 
or not) occurs at lower luminosity at the tip of the RGB 
(L=21308 L and 21252 L for models [1.5M , Z=0.014 
and Z=0.002] respectively). 

- Central helium burning occurs at a higher effective tempera- 
ture, and more extended blue loops are obtained. 

4.1.2. Impact of rotation-induced mixing and thermohaline 
instability 

- As known for a long time (see e.g. lMaeder & Mev net 2000) 
and as shown in Figj2j rotation-induced mixing affects the 
evolution tracks in HR diagram. On the main sequence, ro- 
tational mixing brings fresh hydrogen fuel into the longer 
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lasting convective core, and transports H-burning products 
outwards. This results in more massive helium cores at 
the turnoff than in the standard case and shifts the tracks 
towards higher effec tive temperatures and luminosities all 
along the evolution dEkstrom et aT 20121: Heger & Langerl 



l200fJt iMaeder & Mevnetl 120001: iMevnet & Maederi l2000t) 
The right panel of Figj3]shows the central conditions in stan- 
dard and rotating models for selected masses at solar metal- 
licity. Rotating models behave as star with higher mass dur- 
ing all of their evolution. 

- Thermohaline mixing induced by 3 He-burning becomes 
efficient only on the RGB at the bump luminosit y 
(ICharbonnel & Zahnl [2007bl: ICharbonnel & Lagardell2010l) . 
Beyond this point the double-diffusive instability develops 
in a very thin region located between the hydrogen-burning 
shell and the convective envelope, and has negligible ef- 
fect on the stellar structure. It does not modify the evolution 
tracks in the HR diagram nor in the log T c vs log pc diagram. 

4.2. Lifetimes 

The theoretical lifetimes are shown in Figj4]for the main phases 
of evolution as a function of the initial stellar mass and for the 
four metallicities considered, with the effects of rotation being 
presented. Let us recall the main points: 

- The duration of all the evolutionary phases decreases when 
the initial stella r mass increases and when the rnetal- 
licity decreases (ISchaller et al.1 Il992l ISchaerer et aD 119931: 
ICharbonnel 119941) 7 which is consistent with the above- 
mentioned luminosity and effective-temperature differences. 

- After the bump luminosity the thermohaline mixing does 
not affect the lifetimes on the RGB and on the early-AGB 
phases. 

- The lifetimes of rotating models on the main sequence are 
increased compared to those of the standard models. Indeed 
rotation-induced mixing brings fresh hydrogen fuel in the 
stellar core during that phase. As a consequence, the exhaus- 
tion of hydrogen in the central region is delayed and the life- 
time on the main sequence increases; additionally the mass 
of the He-core is larger at the end of main sequence when ro- 
tation is accounted for (Ekstrom et al.l2012 : Heger & Langerl 



20001: IMaeder & Mevnetll2000l:lMevnet & Maederj|2000l) 



- As a consequence of having a more massive core at the end 
of the main sequence due to rotation, the models that undergo 
the He-Flash spend a shorter time on the RGB (see Fig. @). 
Similarly to the MS, all rotating models have a longer life- 
time during the quiescent central He burning phase. 

- Lifetimes on the early-AGB are longer in rotating models 
than in standard models, due to their more massive core. This 
remains unchanged if we add the lifetime on the TP-AGB 
(1.1. 10 6 yr for [1.5 M , Z ] and 2.10 6 yr for [3 M Q , Z ] ) 
to the early-AGB as the total lifetime on the AGB is only 
increased by a few percent. 

- As seen in previous subsections thermohaline mixing has a 
negligible impact on the stellar structure, and does not mod- 
ify the stellar evolutionary tracks. Similarly it does not im- 
pact on lifetimes. 

5. Global asteroseismic quantities 

In recent years, a large number of asteroseismic observations 
have been obtained for different kinds of stars. In particular, the 
detection and characterization of solar-like oscillations in a large 



numb er of red giants by space missions (e.g. iDe Ridder et al.l 
l2Q09h promises to add valuable and independent constraints 
to current stellar models. The confrontation between models 
including a detailed description of transport processes in stellar 
interiors and these asteroseismic constraints opens a new 
promising path for our understanding of stars. 

Rotation is one of the key processes that change all outputs 
of stellar models (see Sect. [4]i with a significant impact on aster- 
oseismic observables. In the case of main-sequence solar-type 
stars, rotation is found to shift the evolutionary tracks to the blue 
part of the HR diagram resulting in higher values of the large 
frequency separation for rotating m odels than for non-rotatin g 
ones at a given evolutionary stage (lEggenberger et al1l2010al) . 
For red giants, rotating models are found to decrease the 
determined value of the stellar mass of a star located at a given 
luminosity in the HR diagram and to increase the value of 
its age. Consequently, the inclusion of rotation significantly 
changes the fundamental parameters deter mined for a star by 
perform ing an asteroseismic calibration (lEggenberger et all 
[2010bh . 



5.1. Scaling relations and asymptotic quantities provided for 
the present grid 

For this new grid of standard and rotating models, we now pro- 
vide the values of different asteroseismic parameters (see Table 
|2]i that can be directly computed from global stellar proper- 
ties using scaling relations. We also use the information on the 
internal structure of models to compute asymptotic relations. 
These scaling relations are particularly useful to constrain stel- 
lar parameters and to obtain new information about stellar evo- 
lution without the need to pe rform a full asteros eismic analy- 
sis (see e.g. IStello et aljl2009h iMiglio et alj|2009t iHekker et al.l 
2009l:lMosser et al.l20ld:IChaphn et al.l201 lbUBeck et al.l201 1 



Bedding et al.l 1201 ll) . Although solar-like oscillations are ex- 
pected to be excited in relatively cool stars (main-sequence as 
well as red giant stars), these global asteroseismic quantities are 
provided for all models of our grid and not only for models with 
a significant convective envelope. 

The first global asteroseismic quantity we provide is the large 
frequency separation Av, which is expected to be proportion al to 
the square root of the mean stellar density (e.g. lUlrichlll 9861) : 



Av SC2 



scale 



M\ a5 /r- u 



with the solar large frequency separation Av Q = 134.9//Hz 



(8) 



iBelkacem et alJ d201 ll) have proven that the frequency v max 
at which the oscillation modes reach their strongest amplitudes is 
approximatively proportional t o the acous tic cut-off frequency, 
as sus pected by iBrown et al.1 (Il99ll) and iKieldsen & Bedding] 
(fl99l . 



— ^max.0 



MM RY 2 ( r eff v 



with the solar value v maX;0 = 3150^/Hz. 



(9) 



Finally, the value of maximum oscillation amplitude relative 
to that of the Sun (A max /A max0 ) is computed using the following 
relation (e.g. lHuber et al.ll2.Ql ll) : 
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log p c (g cm" 3 ) lQ g Po (g-cm- 3 ) 

Fig. 3. Left panel Central density as a function of central temperature for the rotating Z models. Colours indicate the mass of 
convective regions (convective core and convective envelope) over the total stellar mass. Shaded regions indicate the phases of 
central hydrogen-burning (MS) and helium-burning (He-B). Right panel Central density as a function of central temperature for 
standard and rotating models (dashed and solid lines respectively) for four initial stellar masses from the ZAMS beyond. 




M (M s ) M (M G ) 

Fig. 4. Duration of the different evolution phases as a function of the initial stellar mass, (left) Effects of metallicity for the standard 
models, (right) Effects of rotation at solar metallicity. 



ll 



(M/M Q y\T ( 



- eff 



■ eff.o 



(10) 



with the solar value from Kepler A„ mx B = 2.5ppm. 

Following Huber and collaborators we adopt s=0.838 and 

t=1.32. A value of r = 2 is adopted as advocated by 



iKieldsen & BeddinS d 1995t) . Note that we simply assume 
that these relations hold for the different models computed 
here and we do not test the validity of the scaling relations. 
Preliminary studies suggest that these simple scaling relations 
hold reasonnabry well for main seq uence and RGB stars (e.g. 
IStello et al.ll20(MlWhite et al.ll201 ll) . 
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Fig. 5. Maximal amplitude compared to solar as a function of the large separation: (left panel) at solar metallicity for 1.25 M Q , 
2.0 M , and 6.0 M rotating models, and (right panel ) at two metallicities Z=0.014 (solid lines) and Z=0.0001 (dotted-lines) for 
a 2.0 M rotating model. Evolutionary phases are color-labeled: pre-main sequence (cyan), main sequence (black), RGB (green), 
He-burning (blue), and AGB (red) 



In addition to the asteroseismic observables deduced from 
scaling relations, asymptotic asteroseismic quantities are also 
provided. 



The asymptotic large frequency separation is given by : 



AVasymp 



c s 



(ID 



where R is the stellar radius, and Cj is the sound speed. 



The total acoustic radius T is provided using the following 
relation, which is directly related to the large frequency separa- 
tion : 



Jo c s 2 ■ Av a; 



(12) 



The acoustic radius at the base of convective envelope (Ibce) 
and at the location of helium second-ionisation regiorQ (te e ) are 
determined with following relations : 

r BCE dr r rs ' dr 

fBCE= — ,*He= — , (13) 

Jo C s Jo C s 

where Tbce and rn e represent the stellar radius at the base of con- 
vective envelope and at the location ofhelium second-ionization 
zone. 

The period spacing of gravity modes for l—\ and different 
acoustic radii can be determined with the following asymptotic 
relation : 



7 The Schwarzschild citerium allows us to define the base of the con- 
vective envelope. The mininum in Y\, the adiabatic exponent, corre- 
sponding to He-II defines the location of helium second ionisation re- 
gion. 



AII(£ = 1) = 



(14) 



where N is the Brunt- Vaisala frequency. i\ and X2 define the do- 
main (in radius) where the g modes are trapped. Within this re- 
gion, the mode frequency (a>) must satisfy the following condi- 
tions : 



(J < N 2 



and 



co l <Sf = 



l(l+\)c]. 



(15) 



(16) 



with Si the lamb frequency. 

In the following discussions, we use the large separation de- 
termined with the scaling relation Av sca i e (noted Av in the follow- 
ing sections). We have compared the large separation estimates 
Av S caie (Eq. 8) and Av asymp (Eq. 1 1) to evaluate the difference be- 
tween both expressions. On the main sequence, this difference 
varies between 3 and 5% for 1.0 M model at solar metallicity. 
We find that the relative error obtained when using Eq. |8]or Eq. 
rm depends on the stellar mass (i.e. 3-8% for [2.0 M , Z ] MS 
model), the metallicity (i.e. 10-12% for the [2.0 M , Z=10~ 4 ] 
MS model), and the evolutionnary phase (i.e. <15% on RGB, 
<10% on He-B both for [1.0 M , Z ] model). A forthcoming 
work, out of the scope of the paper, is neces sary to interp r et thes e 
differences and to see how they compare to lWhite et al.l (1201 ll) 

5.2. Evolution of the asteroseismic observables 

We now describe the global asteroseismic properties of low- and 
intermediate-mass stars along their evolution, and discuss the 
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1000 100 



An(i=i) 

Fig. 7. The stellar luminosity as a function of the asymptotic pe- 
riod spacing of g-modes for the standard 2.0 M model at so- 
lar metallicity. Evolutionary phases are color-labeled: subgiant 
(green dashed), red giant (green solid), He-flash episod (orange), 
He-burning (blue), and asymptotic giant branch (red) 



impacts of stellar mass, metallicity, and of non-standard mixing 
mechanisms. We illustrate our discussion with Fig. [5] to [10] In 
Fig|5]we present the maximal oscillation amplitude as a function 
of the large frequency separation all along the evolution from 
the pre-main sequence to the AGB phase for three different ini- 
tial stellar masses at solar metallicity (left panel) and for a 2 M 
star for the lowest (Z=0.0001) and the highest (Z=0.014) metal- 
licities explored here (right panel). The line colors change with 
the stellar evolution phases: cyan corresponds to the pre-main 
sequence, black to the main sequence, green to the red giant 
branch, blue to core helium-burning, and red to the asymptotic 
giant branch. Fig.|6]presents the same quantities for all the stel- 
lar masses considered at solar metallicity, for both the standard 
and the rotating cases (dashed and solid lines respectively), each 
evolutionary phase being presented singly. 

5.2.1 . Trends for a given stellar mass and metallicity - the 
case of the 2.0 M Q , Z models 

In the following we describe the evolution of the seismic proper- 
ties of a star of 2.0 M at solar metallicity as depicted by the solid 
line in Fig|5]and the black line in Fig|6] Although this 2 M Q star 
is not expected to show solar-like oscillations during the main 
sequence due to its too large surface temperature, it represents a 
typical model for a red giant exhibiting solar-like oscillations. 
Moreover, the discussion of the changes of the asteroseismic 
quantities during the main sequence remains globally valid for 
less massive main-sequence stars (the main difference being the 
value of the ratio between the acoustic radius at the base of the 
convective envelope and the total acoustic radius, see Fig. [8]). 

The large separation Av increases along the pre-main se- 
quence because of its inverse dependence in stellar radius, which 



decreases during that phase (see EqJT4li. On the other hand, the 
simultaneous decrease of maximal amplitude A max /A max0 re- 
sults from the decrease of luminosity and from the rise of ef- 
fective temperature (see EqlTOV 

On the main sequence, A max /A max>0 rises because of its de- 
pendence with luminosity and with the inverse of the effec- 
tive temperature. In addition, the expansion of the stellar ra- 
dius causes Av to increase. The radius increase also causes the 
large separation to drop significantly as the model evolves on the 
RGB, while the maximum amplitude increases during this phase 
as a result of the luminosity increase and the effective tempera- 
ture decrease. 

After the He-flash at the tip of the RGB, the luminosity and 
the stellar radius dwindle with a rise of effective temperature, 
implying that the maximum of amplitude wanes while the large 
separation increases. Throughout the helium burning phase, 
the luminosity and stellar radius increase while the effective 
temperature and stellar mass decrease. Consequently, a gain of 
A max is obtained while the large separation reduces. 

Figure [7] presents the asymptotic period spacing of gravity 
modes An(^= 1 ) for the standard mode l of 2.0 M Q at solar metal- 
licity. As proposed by iBedding et al.l d201 ll) and iMosser et al.l 
d201 ll) . this quantity allows to distinguish two stars that have 
the same luminosity, one being at the RGB bump and the other 
one being at the clump undergoing central He burning. Indeed, 
at log(L/L ) ~ 2.0, An(€=l)=55 s at the RGB bump, and 
An(£=l)=190 s in the clump. As the stellar structure, and more 
particularly the presence of the convective core affects the do- 
main where the g-modes a re trapped, An(l= 1 ) is larger i n clump 
stars than in RGB stars (IChristensen-Dals gaardl I20TTI) . Large 
variations of An(£=l) are observed during the thermal pulses 
and He-flash phase because of the formation of an intermediate 
convective zone during the He-flashes and thermal pulses. We 
consider that the modes are trapped in t he outermost radiativ e 
zone leading to a large value of An(^=l) dBildsten et al.ll2012h . 

Figure [8] shows the acoustic radius at the base of convec- 
tive envelope (tBCE) and at the base of He-ionization zone (tn e ), 
both over the total acoustic radius as a function of luminosity, 
for the standard (black line) 1.0 M (left panel) and 2.0 M Q 
(right panel) models at solar metallicity from the main sequence 
to He-burning phase. As the extent (in radius) of the convective 
envelope decreases with increasing stellar mass on the main se- 
quence, the difference between tn e and tecE becomes smaller. 

As the convective envelope deepens inside the star with the 
first dredge-up, the acoustic radius t^cE decreases during the 
subgiant branch. During this phase, the acoustic radius at the 
base of He II ionization zone follows the variation of effective 
temperature and te e increases while T e ff decreases. In addition 
the total acoustic radius increases because of its dependence with 
the stellar radius that increases during this phase. Consequently, 
t He /T decreases from 0.8 to 0.68 and from 0.9 to 0.7 for [1.0 
M Q , Z ] and [2.0 M Q , Z ] models respectively . While the star 
ascends the RGB, tecE follows the convective envelope and the 
acoustic radius tuJT decreases to 0.2 for 2.0 M model until 
the RGB tip is reached (at log L/L ~ 3.5). When the stellar 
luminosity decreases after the RGB tip the convective envelope 
retracts in radius, and the effective temperature increases while 
the star contracts. Therefore tBCE/T and tn e /T increase until the 
central temperature is sufficient to ignite helium. At the end of 
He burning phase, the surface layers expand and the convec- 
tive envelope deepens while the effective temperature decreases. 
During this phase tecE/T and tn e /T in 2.0 M model decrease 
slowly from ~0.62 and ~0.1 to ~0.55 and ~0.05 respectively. 
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Pre-Main sequence Main Sequence 





Fig. 6. Maximal amplitude compared to solar as a function of the larger separation for models at Z . Top-left panel: on the pre-main 
sequence ; top-right panel: on the main sequence ; bottom-left panel: on the red giant branch ; bottom-right panel: on He-burning 
phase. Standard models and rotating models are shown with dashed and solid lines respectively except for pre-main sequence. The 
corresponding evolution in the Hertzsprung-Russell diagram is shown on the right corner of each panel 



Although we do not include this phase in the files presented 
in Sect l3.ll we discuss the evolution of the asteroseismic pa- 
rameters during the TP-AGB. Figure [10] shows the evolution of 
stellar parameters (L,T e ff, and M), and asteroseismic parameters 
(Av, v m a X , and A max ) as a function of time from the first thermal 
pulse for 2.0 M model at Z Q . 



Between each thermal pulse, the stellar radius increases. The 
mass slowly decreases by steps at this phase, remaining almost 
constant during the inter-pulses. Then, the large separation de- 
creases between each thermal pulse as well as v max despite the 
slight decrease of effective temperature. The maximal amplitude 
increases as a result of the luminosity increase and the decrease 
of the effective temperature and total stellar mass. The strong in- 
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Main Sequence Sub-giant RGB He-B Main Sequence Sub-giant RGB He-B 




Log(L/L ) Log(L/L s ) Log(L/L ) log(L/L s ) Log(L/I^) Log(L/L a ) Log(L/I^) log(L/L s ) 



Fig. 8. Acoustic radius at the base of the convective envelope (tBCE, solid line) and at the base of He-ionization zone (te e , dashed 
line) compared to total acoustic radius T as a function of the stellar luminosity for (left panel) 1.0 and (right panel) 2.0 M models 
at solar metallicity following standard (black line) and rotating prescriptions (red line) ; (from left to right) on the main sequence, 
the sub-giant branch, red giant branch, and He-burning phase. 




log Teff (K) log Teff (K) 

Fig. 9. Color-coded HR diagram for all the masses in the grid at solar metallicity. The color code represents the values (in seconds) 
of the acoustic radius at the base of the ionization of He II fee (left panel) and of the asymptotic period spacing of g-modes AYl(£ = 1 ) 
(right panel). The values increase from blue to red in both cases as shown on the scales on the right of the plots. 

crease of A max during the interpulses is due to the large mass loss 5.2.2. Effects of metallicity 
at this phase, better seen in the last pulse shown in Fig. [10] 

Although the current asteroseismic missions focus on solar 
metallicity stars, we present the effects of metallicity on the 



12 



N. Lagarde et al.: III. Grid of stellar models and asteroseismic quantities 



global asteroseismic parameters; they are depicted in Figj5](right 
panel) for a 2.0 M Q model. 

As discussed in l4.1.1l at lower metallicity the main sequence 
is shifted to the blue and to higher luminosity on HR diagram. 
Therefore, the track in the A max vs Av plot is moved to lower 
A max /A m£IV0 values. In addition, the luminosity of He-ignition at 
the RGB tip occurs at lower A max /A mflt0 when the metallicity is 
lower. The track all along the evolution in the A max vs Av plot is 
shifted to the right (toward higher Av values) when the metallic- 
ity decreases. This is a direct consequence of a more metal-poor 
star behaving as a more massive star. 

5.2.3. Effects of stellar mass 

When the mass increases, the ZAMS (point at higher Av and 
lower A max ) is shifted to lower Av and higher A max (see Fig. [5] 
and top panels of Fig|6]l. Indeed, the large frequency separation 
at the point of ZAMS decreases because of its dependence with 
the mean stellar density. Similarly to the Hayashi tracks during 
the pre-main sequence, stars with different masses tend to join 
together in the A max vs Av plot during the RGB and early-AGB 
phase (see green and red lines respectively in Fig. |5). When 
stellar mass increases stars follow more extended blue loops 
during the central He-burning phase. As shown in Fig.|6](bottom 
right panel), blue loops range over high Av and low A max /A max>0 . 

Figure|9]shows the theoretical evolution of the stellar acous- 
tic radius at the base of He-ionization zone and of the asymptotic 
period spacing of g-modes along the evolutionary tracks in the 
HR diagram for models of different masses at solar metallicity. 
Both quantities variations are expressed in seconds. Evolution 
is shown from the pre-main sequence to the early-AGB phase 
and from the subgiant branch to the early-AGB phase for tn e 
and An respectively. At a given evolutionary phase, the larger 
the stellar mass, the larger also the acoustic radius at the base of 
He-ionization zone. The same is true for An. 



5.2.4. Impact of rotation on the evolution along the amplitude 
vs Av diagram 

In this section, we focus on the effect of accounting for rotational 
mixing on the global asteroseismic parameters. As thermoha- 
line mixing solely affects the abundance at the surface and at 
the external layers of HBS pattern during the giant phase, with- 
out modifying the temperature, radius, luminosity or mass of the 
models, this process has no impact on the global asteroseismic 
parameters. It is nevertheless not excluded that well chosen as- 
teroseismic parameters might help constraining the thermoha- 
line mixing. 

As discussed in Sect. [4] rotation modifies the position of the 
evolution tracks in the HR diagram (see Fig. Consequently, 
differences appear on the maximal amplitude and the large sepa- 
ration, between standard and rotating models. This effect clearly 
shows up at the end of the main sequence due to the increased 
width of the main sequence in rotating models. 
Then on the subgiant branch rotating models evolve at higher lu- 
minosities than their standard counterparts, which implies higher 
maximal amplitude. Rotation has no impact on the ratios tBCE/T 
and tne/T as seen in Fig. [8] In this figure, the effect of rotation on 
the stellar luminosity during the subgiant phase is also shown. 
When the stars reach the RGB, the differences in effective tem- 
perature and luminosity between standard and rotating models 
become marginal; therefore the tracks are almost identical in the 
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Fig. 10. From top to bottom: Evolution of global stellar 
properties (luminosity L, effective temperature T e ff and to- 
tal stellar mass M) and asteroseismic parameters (Av, and 
v max ,A max /A maXj0 ) on the TP-AGB of the [2.0 M , Z ] model 
computed with rotation up to the RGB tip and thermohaline mix- 
ing all along the evolution. The abscissa is the time since the first 
thermal pulse. 



Figure [6] (bottom left panel). 

Rotation affects the extension and width of the blue loops. As 
a consequence, the differences in asteroseismic parameters ap- 
pear during the combustion phase of He, as illustrated in Fig. [6] 
(bottom right panel). 

5.3. Same point on HR-diagram, same asteroseismic 
parameters ? 

Figure QT| shows how asteroseismic parameters can help to dis- 
tinguish two stars with the same effective temperature and lu- 
minosity. On the left panel we show the 1.5 M and 2.0 M Q 
tracks in the HR diagram with the intersection indicated by the 
black dot. At this point, the 1.5 M star is evolving on the HR 
gap and has the same effective temperature and luminosity as 
the 2.0 M Q star evolving on the pre-main sequence. This black 
dot corresponds to two different points in the A max vs Av plot 
shown on the right panel. As we can see, the larger separation of 
the 1.5 M model is lower that of the 2.0 M , while the maxi- 
mal amplitude is higher. These stars do not have the same mean 
stellar density, and therefore the same Av. On the other hand, the 
total stellar mass is different at this point, and hence the maximal 
amplitude. 



6. Conclusions 

In this study we present a grid of single star evolution 
models in the mass range between 0.85 M e and 6.0 M , 
for four metallicities, including the impact of rotation 
and thermohaline mixing, along with standard models. 
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log Teff (K) Au (//.Hz) 

Fig. 11. Left panel : Evolution tracks in Hertzsprang-Russel dia- 
gram of 1.5 and 2.0 Z (solid and dotted lines respectively) stars 
at solar metallicity. Right panel : Maximal amplitude over sun 
value as a function of the large separation for the same models. 
As in figure [5] the lines color indicates the evolutionary phase 

All data detailed in Table [2] are available on the website 
http://obswww.unige . ch/Recherche/evol/-Database- 
for all the models computed (i.e., standard on one hand, and 
with rotation-induced mixing and thermohaline instability on 
the other hand) . 

We recall the various impacts of metallicity variations and 
rotation-induced mixing on the stellar properties that were 
already discussed in literature. Thermohaline mixing does not 
change the stellar parameters like luminosity and effective 
temperature. As such this process does not affect the seismic 
properties as analyzed here. However and as discussed in Paper 
I and II, it changes the surface abundances from the bump 
luminosity on, and has to be taken into account to explain the 
observed chemical properties of bright giant stars. 

Last but not least we also present the evolution of global 
asteroseismic parameters for all the models in our grid. The 
large frequency separation Av, the frequency v max , and the maxi- 
mum oscillation amplitude A max , are computed using scaling re- 
lations. Asymptotic asteroseismic quantities are also computed 
: Avasympt, another estimate or the large frequency separation; 
?bce, the acoustic radius at the base of convective envelope; tan, 
the acoustic radius at the base of the Hell ionization zone ; T, 
the total acoustic radius; and An(^=l) , the period spacing of 
gravity modes. 

We show that rotation-induced mixing has an impact on 
these quantities, contrary to thermohaline mixing. While rota- 
tion changes the global properties of main sequence stars and 
has an impact on the global asteroseismic properties, thermoha- 
line mixing is negligible on these aspects although it changes the 
surface abundances in the red giants. In addition to spectropho- 
tometric studies, seismic studies allow to distinguish two stars 



with approximatively the same luminosity and effective temper- 
ature but with different evolutionary stages. 
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